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ABSTRACT 



Context. DBV stars are pulsating white dwarfs with atmospheres rich in He. Asteroseismology of DBV stars can provide valuable 
clues about the origin, structure and evolution of hydrogen-deficient white dwarfs, and may allow to study neutrino and axion physics. 
Recently, a new DBV star, KIC 8626021, has been discovered in the field of the Kepler spacecraft. It is expected that further monitoring 
of this star in the next years will enable astronomers to determine its detailed asteroseismic profile. 

Aims. We perform an asteroseismological analysis of KIC 8626021 on the basis of fully evolutionary DB white-dwarf models. 
Methods. We employ a complete set of evolutionary DB white-dwarf structures covering a wide range of effective temperatures and 
stellar masses. They have been obtained on the basis of a complete treatment of the evolutionary history of progenitors stars. We 
compute g-mode adiabatic pulsation periods for this set of models and compare them with the pulsation properties exhibited by KIC 
8626021. 

Results. On the basis of the mean period spacing of the star, we found that the stellar mass should be substantially larger than 
spectroscopy indicates. From period-to-period fits we found an asteroseismological model characterized by an effective temperature 
much higher than the spectroscopic estimate. 

Conclusions. In agreement with a recent asteroseismological analysis of this star by other authors, we conclude that KIC 8626021 is 
located near the blue edge of the DBV instability strip, contrarily to spectroscopic predictions. We also conclude that the mass of KIC 
8626021 should be substantially larger than thought. 

Key words, stars: evolution — stars: interiors — stars: white dwarfs 
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1. Introduction 

V777 Her (or DBV) stars are ^-mode variable white dwarfs 
with He-rich atmospheres (DB) and effective temperatures in 
the range 21 500 < Teff < 29 000 K that pulsate with periods 
between 100 and 1100 s (Winget & Kepler 2008; Althaus et al. 
2010). They are the hotter cousins of the ZZ Ceti (or DAY) stars, 
which are pulsating H-rich atmosphere white dwarfs that define 
an instability strip centered at T^ff * 12 000 K. As an excellent 
demonstration of the validity of the stellar pulsation theory, the 
existence of the DBV class of compact pulsators was predicted 
on theoretical grounds (Winget et al. 1982a) before it were con- 
firmed observationally (Winget et al. 1982b). Pulsations in V777 
Her are thought to be driven by a combination of the K-y mech- 
anism acting in the He partial ionization zone — and thus setting 
the blue edge of the DBV instability strip (Winget et al. 1983; 
Bradley & Winget 1994; Corsico et al. 2009), and the "convec- 
tive driving" mechanism (Brickhill 1991; Goldreich & Wu 1999) 
which is thought to be dominant once the outer convection zone 
has deepened enough. 

White-dwarf asteroseismology — the comparison between 
the pulsation periods of white dwarfs and the periods computed 
for appropriate theoretical models — allows us to infer details of 
the origin, internal structure and evolution of white dwarfs. In 
particular, estimates of the stellar mass. He and H layer mass, 
core composition, magnetic field, rotation rate, seismological 
distance, and properties of the outer convection zone can be 



inferred from the observed pulsations of DAV and DBV stars 
(Winget & Kepler 2008; Althaus et al. 2010). Finally, an even- 
tual measurement of the temporal changes in the observed stable 
periods of DB Vs could allow to study neutrino emission (Winget 
et al. 2004) and also axion emission (see Isern et al. 2010 and 
references therein). 

Unlike the ZZ Ceti variables, that constitute the most nu- 
merous class of pulsating white dwarfs (about 150 objects are 
known today; Castanheira et al. 2010), V777 Her are quite rare 
and difficult to find, and until recently, only 20 stars of this class 
were known (Beauchamp et al. 1999; Nitta et al. 2009; Kilkenny 
et al. 2009). The list is now a bit enlarged with KIC 8626021 
(= GALEX J192904.6-H444708, T^ff = 24950 + 750 K and 
log ^ = 7.91+0.07 dex), a very recently discovered DBV star lo- 
cated in the field of view of the Kepler Mission by 0stensen et al. 
(201 1) (hereinafter 0EA1 ifl This star exhibits at least five peri- 
odicities in the range 197 - 376 s and amplitudes between 1 and 
5.2 mma, with the three strongest modes showing a triplet struc- 
ture due possibly to rotation. It is expected that additional modes 
with higher radial orders (longer periods) will be detected in fu- 
ture runs of Kepler, thus increasing the potential of asteroseis- 
mology to infer its interior structure. Given the values of T^s and 
log g quoted above and the DB models of Althaus et al. (2009a), 
KIC 8626021 should have a mass of M. = 0.56 ± O.O3M0. 



' Even more recently, Hermes et al. (201 1) have reported the discov- 
ery of the first DAV star in the Kepler Mission field. 
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Motivated by the exciting discovery of the first pulsating 
white dwarf in the Kepler field of view, we present an astero- 
seismological analysis of KIC 8626021 on the basis of the DB 
white-dwarf models presented in Althaus et al. (2009a). This 
grid of evolutionary models was computed for a wide range of 
stellar masses on the basis of a complete treatment of the evolu- 
tionary history of progenitors stars, including those stages rele- 
vant for the chemical profiles of the white dwarf, such as the core 
H and He burning phases, the thermally pulsing asymptotic giant 
branch phase, and the born-again episode that is responsible for 
the H deficiency. 

By considering the mean period spacing exhibited by the 
star, we found that KIC 8626021 should be substantially more 
massive than spectroscopy indicates. We also found that period 
to period fits favour an asteroseismological solution character- 
ized by an effective temperature much higher than the spectro- 
scopic estimate. 

While writing this paper, BischofF-Kim & 0stensen (2011) 
(hereinafter BK011) announced the results of their own inde- 
pendent asteroseismological analysis on KIC 8626021. They 
found that this star is actually a hot DBV. Using entirely different 
models and methods, we arrive at the same conclusion. The fun- 
damental difference between our approach and that of BK01 1 is 
that we run evolutionary models starting on the Zero Age Main 
Sequence, accounting for nuclear burning, time dependent diffu- 
sion of the elements and other physical processes that take place 
in the course of the evolution of a star BK011 perform "fast 
white dwarf asteroseismology", where they guess and parame- 
terize the internal chemical composition profiles and build static 
models. That method allows them a fuller exploration of param- 
eter space. For KIC 8626021, they find very good fits (hotter 
than what the preliminary spectroscopy suggests) but also con- 
clude that their internal Oxygen composition profiles are in poor 
agreement with stellar evolution calculations. They suggest fur- 
ther studies to see what the physical parameters of models that 
agree with stellar evolution calculations would be. In essence, 
this is the kind of study we present here. While our best fit mod- 
els are evidently different, we concur with the conclusion that 
KIC 8626021 is hotter (and more massive) than suggested by 
the initial spectroscopic study. The conclusion that this star is 
residing at the blue edge of the DBV instability strip appears 
to be robust, and calls for the necessity of a new and improved 
spectroscopic determination of its effective temperature. 

In Section |2] we present some details about our models and 
methods. In Section[3] we use the average period spacing of KIC 
862602 1's period spectrum to draw some conclusions about its 
mass and effective temperature. We contrast our models to the 
grid of models computed by BK011. In Section H) we present 
our best fit models. Again, we contrast our results with those of 
BK011. In Section |5] we discuss the discrepancies between the 
results coming from spectroscopy and from asteroseismology, 
and the differences between our asteroseismological results and 
those of BK01 1. We conclude in Section |6] 

2. Modeling 

2. 1 . Numerical codes and asteroseismological approach 

The evolutionary DB white dwarf models employed in this work 
were presented in Althaus et al. (2009a) and we refer to that 
paper for details. Briefly, the models were computed with the 
LPCODE stellar evolutionary code we employed in our pre- 
vious studies on the formation of PG 1159 and extreme hor- 
izontal branch stars (Althaus et al. 2005; Miller Bertolami & 
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Fig. 1. The location of the known DBV stars on the T^g - \ogg 
plane. Also included are our DB white dwarf evolutionary tracks 
displayed with different colors according to the stellar mass. The 
location of KIC 8626021 according to spectroscopy is empha- 
sized with a red star symbol, while the location of our astero- 
seismological models for this star is marked with orange circles 
connected with a dotted line. The largest circle corresponds to 
our best-fit model (model 2; see Sect. HJl. For comparison, we 
include the location of the best-fit asteroseismological model of 
BK011 with a blue star symbol. The theoretical blue edge of 
the DBV instability strip corresponding to the version MLT2 
(a - 1.25) of the MLT theory of convection, as derived by 
Corsico et al. (2009), is depicted with a blue dashed line. Note: 
the DBV stai- EC 05221 (Kilkenny et al. 2009) is not included 
because at present there is no any estimation of T^s and log g 
available. 



Althaus 2006; Miller Bertolami et al. 2008), hot DQ white 
dwarfs (Althaus et al. 2009b) as well as the formation and evo- 
lution of He-core white dwarfs with high metallicity progenitors 
(Althaus et al. 2009c). The LPCODE evolutionary code consid- 
ers a simultaneous treatment of non-instantaneous mixing and 
burning of elements, which is of primary importance for the cal- 
culation of chemical abundance changes during the short-lived 
evolutionary stages characteristic of unstable burning episodes, 
like the born-again stage, from which our starting H-deficient 
white dwarf configurations are derived. Nuclear reaction rates 
are from Caughlan & Fowler (1988) and Angulo et al. (1999). 
The '^C(q', y)"'0 reaction rate was taken from Angulo et al. 
(1999), which is about twice as large as that of Caughlan & 
Fowler (1988). A moderate, diffusive overshooting in the core 
and in the envelope is allowed during pre-white-dwarf evolution. 
For the white dwarf regime, we considered the following main 
physical ingredients. Neutrino emission rates for pair, photo, 
and bremsstrahlung processes are those of Itoh et al. (1996). 
For plasma processes, we use the treatment presented in Haft 
et al. (1994). Radiative opacities are those of the OPAL project 
(Iglesias & Rogers 1996), including C- and 0-rich composition. 
We adopted the conductive opacities of Cassisi et al. (2007). This 
prescription covers the whole regime where electron conduction 
is relevant. For the high density regime, we used the equation of 
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state of Segretain et al. (1994), which accounts for all the impor- 
tant contributions for both the liquid and solid phases. For the 
low-density regime, we used an updated version of the equation 
of state of Magni & Mazzitelli (1979). Convection was treated in 
the formalism of the mixing-length theory as given by the ML2 
parameterization (Tassoul et al. 1990). All of our white dwarf 
sequences were computed in a consistent way with the evolution 
of the chemical abundance distribution caused by element diffu- 
sion along the whole cooling phase. In particular, we considered 
gravitational settling and chemical diffusion of "^He, '^C, '■'C, 
^^N, and '^O. 

The adiabatic pulsation periods employed in the present 
study were assessed with the help of the pulsational code de- 
scribed in Corsico & Althaus (2006). The prescription we fol- 
low to assess the run of the Brunt-Vaisala frequency (A^) is 
the so-called "Ledoux Modified" treatment — see Tassoul et al. 
(1990) — appropriately generalized to include the effects of hav- 
ing three nuclear species (O, C, and He) varying in abundance. 

Our asteroseismological approach basically consists in the 
employment of detailed white dwarf models characterized by 
very accurate physical ingredients. These models are obtained 
by computing the complete evolution of the progenitor stars. 
We have applied successfully this approach to the hot DOVs 
or GW Vir stars (see Corsico et al. 2007a, 2007b, 2008, 2009) 
and recently to an ensemble of bright ZZ Ceti stars (Romero et 
al. 2011). Since the final chemical stratification of white dwarfs 
is fixed in prior stages of their evolution, the evolutionary his- 
tory of progenitor stars is of utmost importance in the context 
of white dwarf asteroseismology. Our asteroseismological ap- 
proach, while being physically sound, is by far more computa- 
tionally demanding than other approaches in which simplified 
models are used. As a result, our approach severely limits the 
exploration of the parameter space of the models. Indeed, for the 
case of DB white dwarfs, we have only two parameters which 
we are able to vary in a consistent way: the stellar mass (M») and 
the effective temperature (Teff). Instead, the content of He (Mue), 
the shape of the C-0 chemical structure at the core (including the 
precise proportions of central O and C), and the thickness of the 
chemical transition regions are fixed by the evolutionary history 
of progenitor stars. 

2.2. DB white dwarf evolutionary models 

The initial models for our DB white dwarf sequences corre- 
spond to realistic PG 1159 stellar configurations derived from 
the full evolutionary calculations of their progenitor stars (Miller 
Bertolami & Althaus 2006). All the sequences were computed 
from the ZAMS through the thermally-pulsating and mass-loss 
phases on the AGB and finally to the born-again stage where the 
remaining H is violently burnt. After the born again episode, the 
H-deficient, quiescent He-burning remnants evolve at constant 
luminosity to the domain of PGl 159 stars with a surface chemi- 
cal composition rich in He, C and O (Miller Bertolami & Althaus 
2006). This new generation of PG 1 159 evolutionary models has 
succeed in explaining both the spread in surface chemical com- 
position observed in most PG 1 159 stars and the location of the 
GW Vir instability strip in the log Teff - log ^ plane (Corsico et al. 
2006). Also, these PG 1159 models have been employed in de- 
tailed asteroseismological studies of six pulsating PGl 159 stars 
(Corsico et al. 2007a,b, Corsico et al. 2008, 2009). 

Specifically, we considered nine DB white 
dwarf sequences with stellar masses of: 
0.515, 0.530, 0.542, 0.565, 0.584, 0.609, 0.664, 0.741, and 
O.87OM0. These DB sequences are characterized by the max- 
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Fig. 2. The internal chemical structure (upper panel), and the 
squared Brunt-Vaisala and Lamb frequencies for ( - \ (lower 
panel) corresponding to our template DB white dwarf model 
with a stellar mass - O.565M0, an effective temperature 
Teif ~ 28 400 K, and a He envelope mass of Mhs/M. ~ 6.7x10"^ 
(thick lines). For comparison, we show with thin lines the chemi- 
cal profiles and the Brunt-Vaisala frequency of the best fit model 
for KIC 8626021 found by BK01 1. 



imum He-rich envelope that can be left by prior evolution if 
we assume that they are the result of a born-again episode. 
The value of envelope mass ranges from Mue/Mt ~ 2 x 10"^ 
(M, = 0.515Mo) to MnJM, ~ 1 x 10"^ (M, = O.87OM0). 
The complete set of DB white dwarf evolutionary sequences is 
displayed in Fig.[T]with different colors according to the value of 
the stellar mass. Also shown is the location of the known DBV 
stars with values of Teff and log g extracted from Beauchamp et 
al. (1999) (for the case in which no traces of H are considered; 
black circles), from Nitta et al. (2009) for the nine DBV stars 
of the SDSS (black squares), and Kilkenny et al. (2009) (black 
triangle) . The location of KIC 8626021 according to 0EA1 1 is 
highlighted with a red star symbol. 

2.3. A template model 

We briefly describe the main pulsation properties of our DB 
white dwarf models. We focus on a template DB model with 
M, = 0.565Mo and T^ff ~ 28 400 K. In Fig.|2]we depict the in- 
ternal chemical structure of such model (upper panel), where the 
abundance by mass of the main constituents ("^He, '^C, and '^O) 
is shown in terms of the outer mass fraction [- log(l - Mr/M,)]. 
The chemical structure of our models consists of a C/O core — 
resulting from the core He burning of the previous evolution — 
shaped by processes of extra mixing like overshooting. The core 
is surrounded by a mantle rich in He, C, and O which is the 
remnant of the regions altered by the nucleosynthesis during 
the thermally pulsing asymptotic giant branch. Above this shell, 
there is a pure He mantle with a mass Mue/Mt ~ 6.7 x lO"-', con- 
structed by the action of gravitational settling that causes He to 
float to the surface and heavier species to sink. A double-layered 
structure of the He-rich envelope is clearly visible (Dehner & 
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Kawaler 1995; Gautschy & Althaus 2002; Althaus & Corsico 
2004). In our models, the shape of the double-layered struc- 
ture evolves with time (with appreciable changes even within 
the DBV instability strip) by virtue of time-dependent chemical 
diffusion processes (see Althaus & Corsico 2004). 

The lower panel of Fig. |2] displays the run of the two criti- 
cal frequencies of nonradial stellar pulsations, that is, the Brunt- 
Vaisala frequency and the Lamb frequency (Lf) for { - 1. 
The precise shape of the Brunt-Vaisala frequency largely de- 
termines the properties of the ^-mode period spectrum of the 
model. In particular, each chemical gradient in the model con- 
tributes locally to the value of A^. The most notable features 
are the very peaked feature at the C/O chemical transition 
[-log(l - M,./M,) ~ 0.3], and the bump at the He/C/O inter- 
face [-log(l - M,./M,) ~ 1.5]. On the other hand, the He/C/O 
transition region at -log(l - M, /M») ~ 5 is very smooth and 
does not affect the pulsation spectrum much. For comparison, we 
also include in Fig. |2]the core chemical profiles and the Brunt- 
Vaisala frequency of the best fit model for KIC 8626021 found 
by BK011 (thin lines). 

3. Period spacing: the global information 

For ^-modes with high radial order k (long periods), the separa- 
tion of consecutive periods i\Ak\ - I) becomes nearly constant at 
a value given by the asymptotic theory of nonradial stellar pul- 
sations. Specifically, the asymptotic period spacing (Tassoul et 
al. 1990) is given by: 



An^ = Do/ ^fi{FTTj, 



where 



Jr, 



r 



(1) 



(2) 



This expression is rigorously valid for chemically homogeneous 
stars. In principle, one can compare the asymptotic period spac- 
ing computed from a grid of models with different masses and 
effective temperatures with the mean period spacing exhibited by 
the star, and then infer the value of the stellar mass. This method 
has been applied in numerous studies of pulsating PG 11 59 stars 
(see, for instance, Corsico et al. 2007a,b, 2008, 2009 and ref- 
erences therein). For the method to be valid, the periods exhib- 
ited by the pulsating star must be associated with high order g- 
modes, that is, the star must be within the asymptotic regime of 
pulsations. KIC 8626021 exhibits pulsation periods associated 
to low order modes, with presumably k < 10. Furthermore, the 
interior of DB white dwarf stars are supposed to be chemically 
stratified and characterized by strong chemical gradients built up 
during the progenitor star life. So, the direct application of the 
asymptotic period spacing to infer the effective temperature and 
stellar mass of KIC 8626021 may appear somewhat question- 
able. However, the observed period spectrum is roughly evenly 
spaced, indicating that using asymptotic theory may not be an 
unreasonable thing to do. In order to compare our models to the 
period spacing of the observed pulsation spectrum, we calculate 
the average of the computed period spacings, using: 



Anthe(M., Jeff) = 7 T- y AHi, 

(n - 1) ^ 



(3) 




Q 

5 



i: 



0.56 



24000 25000 



I 



where AII^ is the "forward" period spacing defined as AII^ - 
Ylk+i - Iljt, and n is the number of computed periods laying in the 



Fig. 3. The location of the models of BK011 with an aver- 
age period spacing close to that of KIC 8626021 in the mass- 
temperature plane. The colour scale is related to the the differ- 
ence between the observed average period spacing of the star 
and the period spacing computed with Eq. Q. The spectroscopic 
values for the mass and temperature of the star are also indicated 
with the point and error bars. It is obviously outside the region 
of the best fits in average period spacing. The white regions are 
regions where models failed to converge. They do not affect the 
results. 



range of the observed periods. The theoretical period spacing of 
the models as computed through Eq. ([TJ and (O share the same 
general trends (that is, the same dependence on and Teff), al- 
though An^ is usually somewhat higher than Allthe provided that 
this last quantity is computed on a range of periods associated to 
low order modes (short periods). We estimate the observed mean 
period spacing of KIC 8626021 by means of a nonlinear least- 
squares fit by considering the five periods observed in the star 
(first column of Table|2]i. We obtain AlTobs = 35.78 + 0.47 s. 

To check the validity of using the period spacing to study 
KIC 8626021, we use the grid of models computed by BK011 
(400,000 models). For the fixed internal chemical composition 
profiles depicted in Fig. |2] and for each model, we calculate the 
difference between the observed average period spacing of the 
star and the calculated period spacing [Eq. (O]. The result is 
presented in Fig. [3] A striking feature is the diagonal region 
where the spacings are close. This is a direct consequence of 
Eq. (|2]i and evidence that asymptotic theory applies at least in 
broad strokes. In Eq. (|2]i, the dependence on the Brunt-Vaisala 
frequency is such that the asymptotic period spacing is greater 
when the mass and/or temperature of the model is higher. This 
dependence is apparent in Fig. [3] In order to keep a constant av- 
erage period spacing (that of the star), one has to compensate 
lower temperatures with higher masses and vice versa. 

In Fig. ID we show the run of average of the computed pe- 
riod spacings [Eq. Q] in terms of the effective temperature for 
all of our DB white-dwarf evolutionary sequences. The curves 
displayed in the plot are somewhat jagged. This is because the 
average of the computed period spacings is evaluated for a fixed 
period interval, and not for a fixed A:-interval. As the star evolves 
towards lower effective temperatures, the periods generally in- 
crease with time. At a given Tetf, there are « computed periods 
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Fig. 4. The average of the computed period spacings corre- 
sponding to our DB white dwarf sequences with different stel- 
lar masses. The location of KIC 862 602 1 is shown with a red 
star symbol {T^ff = 24900 + 750 K, ADobs = 35.78 + 0.47 s). 
By simple linear interpolation, we found that the mass of the 
star according to its period spacing is M, = 0.696 ± 0.03 IMq. 
According to spectroscopy (logg = 7.91 + 0.07; see Fig.[TJ, the 
period spacing of the star should be somewhere within the dark 
shaded area. 



laying in the chosen period interval. Later, when the model has 
cooled enough, it is possible that the accumulated period drift 
nearly matches the period separation between adjacent modes 
(|AA:| - 1). In these circumstances, the number of periods laying 
in the chosen (fixed) period interval is n + 1, and Allthe exhibits 
a little jump. In order to iron out the curves of Allthe, in con- 
structing Fig.|4]we have considered pulsation periods in a wider 
range of periods (100 < lit < 1200 s) than that shown by KIC 
8626021 (195 < < 380 s)B The location of KIC 8626021 is 
emphasized with a red star symbol. If we adopt the spectroscopic 
effective temperature of KIC 8626021 (reft = 24900 + 750 K; 
0EA11) we infer a mass for the star according to its mean pe- 
riod spacing of M» = 0.696 + 0.03 IMq, which is ~ 24% larger 
than the spectroscopic estimate (M* = 0.56 + 0.03Mq). The 
dark shaded area in the figure shows the location of the star 
(Allthe, 7"eff) according to the constraint given on gravity by spec- 
troscopy (log^ = 7.91+0.07; see Fig.[TJ. The strong discrepancy 
between the predictions of spectroscopy and the period spacing 
estimator is vividly depicted in this figure. This difference be- 
tween the spectroscopic and seismic masses is particularly rel- 
evant as, ultimately, both estimates of the stellar mass of KIC 
8626021 are based on the same DB white dwarf evolutionary 
models. 

This estimate of the stellar mass is based on the reliability 
of the spectroscopic determination of the effective temperature. 
In this sense, it would be quite interesting to have an estimate 



^ If we adopt a shorter range of periods, closer to the range of periods 
exhibited by KIC 8626021 (say 195 - 380 s), the curves we obtain are 
much more irregular and jumped, although the results in relation to KIC 
8626021 do not appreciable change. 



of the stellar mass devoid of any possible uncertainty from the 
spectroscopic analysis. If we relax for the moment the constraint 
imposed by the effective temperature, and look for the possible 
range of masses for the star to be within the DBV instability strip 
(21 500 < Teff < 29 000 K), we find that 0.60 < MJMq < 0.87. 
Confirming the findings of BK01 1 (see also Fig. |3]l, the rather 
low mass value predicted by spectroscopy (M, ~ O.56M0) is 
excluded by the mean period spacing, unless the star is outside 
the instability strip (Teff > 32 000 K; see Fig. |4]i. This conclu- 
sion relies only on the measured mean period spacing, which is 
a quantity accurately measured, and on the average of the com- 
puted period spacing, a global quantity that is not affected by the 
precise shape of the chemical transition regions of the models. 
It is particularly worth noting that this result is independent of 
the mass of the He-rich envelope. Indeed, Tassoul et al. (1990) 
have shown that the asymptotic period spacing of DBV white 
dwarfs is particularly insensitive to the mass of the He-rich man- 
tle, changing less than 1 s for a change of 10 orders of magnitude 
in the mass of the He-rich mantle (see their Figure 42). Also, this 
conclusion is independent of the effective temperature estima- 
tion of the star So, the result that KIC 8626021 is more massive 
and/or hotter than indicated by spectroscopy seems to be a robust 
finding 



4. Period fits: letting the pulsation periods to speak 
by themselves 

Another way to infer the stellar mass, effective temperature, and 
details of the internal structure of KIC 8626021 is through its 
individual pulsation periods, which are extracted from 0EA11 
and are listed in the first column of Table |2] In this approach we 
seek a pulsation DB white dwarf model that best matches the 
pulsation periods of KIC 862602 1 . We assume that all of the ob- 
served periods correspond to ^' = 1 modes, to be consistent with 
the mean period spacing. The goodness of the match between 
the theoretical pulsation periods (11^^) and the observed individ- 
ual periods (Ilobs,;) is measured by means of a quality function 
defined as: 



1 

-y 



min[(n, 



obs,(' 



n,)2], 



(4) 



where (- 5) is the number of observed periods. The 
DB white dwarf model that shows the lowest value of 
is adopted as the "best-fit model" (see Corsico et 
al. 2007a,b, Corsico et al. 2008, 2009, Romero et al. 
2011). We evaluate the function x^iMt,Tf.ff) for stellar 
masses of 0.515,0.530,0.542,0.565,0.584,0.609,0.664,0.741, 
and 0.870Mq. For the effective temperature we employed a much 
more finer grid (Areff = 10-30 K). The quality of our period 
fits is assessed by means of the average of the absolute period 
differences, 6 - (E/Ii where 6, - Hobs,; - 11^, and by the 

root-mean-square residual, cr = V(2 l^/P)/-^ - V?- 

The quantity (x^y^ in terms of the effective temperature for 
different stellar masses is shown in Fig.|5]together with the spec- 
troscopic effective temperature of KIC 8626021 (red line) and 
its uncertainties (gray strip). We found one strong maximum for 
a model with M. = O.664M0 and Tes = 27 260 K (model 2). 



^ If we apply the method described above but employing the asymp- 
totic period spacing Allj! [Eq. (T)], we obtain M, ~ O.SOM©, and if 
we relax the constraint of the spectroscopic T^ff, we find that 0.66 < 
MJMq < 0.90 for the star to be within the DBV instability strip. 
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Table 1. Asteroseismological solutions for KIC 8626021. 

















model 




^ogg 


M, 


6 




BIC 




[K] 


[cm s"-] 


[Me] 


[s] 


[s] 




1 


29441 


7.993 


0.609 


2.314 


2.721 


1.15 


2 


27 263 


8.099 


0.664 


1.582 


1.934 


0.85 


3 


24 856 


8.235 


0.741 


3.334 


3.685 


1.41 



Such a pronounced maximum in the inverse of implies an 
excellent agreement between the theoretical and observed peri- 
ods. Notably, the effective temperature of this model is much 
higher (about 2400 K larger) than the spectroscopic effective 
temperature of KIC 8626021. Another maximum, albeit much 
less pronounced, is encountered for quite hotter and somewhat 
less massive model with Teff - 29440 and M, = O.6O9M0 
(model 1). According to its Tetj, this model is in excellent agree- 
ment with the asteroseismological solution found by BK011 
(Teff = 29 200 K; see Table [3] below). Finally, a third astero- 
seismological solution is found at a model with M, = O.741M0 
and reft = 24 856 (model 3). Interestingly, this model has ex- 
actly the spectroscopically inferred effective temperature of KIC 
8626021, Teft ^ 24900 K. In Table [1] we summarize the main 
characteristics of the three solutions we found in our analysis, 
that is, models 1, 2 and 3. Models 1 and 3 constitute acceptable 
asteroseismological solutions. However, because the agreement 
between observed and theoretical periods for these models are 
much poorer than for model 2 (see columns 5 and 6 of Table 
[TJ, we adopt this last model as the best-fit asteroseismological 
model of KIC 8626021. One should bear in mind that the anal- 
ysis of the mean period spacing suggests a preferred mass of 
~ O.7M0, between our 0.664 and 0.741Mo models. Thus a 
better fit, with a mass between 0.664 Mq and O.741M0, might 
had arisen if a thinner mass grid were available. 

A detailed comparison of the observed m-Q periods in KIC 
8626021 with the theoretical periods of the best-fit asteroseismo- 
logical model is provided in Table |2] For this model, we obtain 
6 - 1.582 s and <j - 1.934 s. The mean period spacing of our 
best fit model is AO = 35.10 + 0.50 s (nonlinear least-squares 
fit), in excellent agreement with the mean period spacing of KIC 
8626021 (Allobs = 35.78 + 0.47). 

BK01 1 perform detailed period fits to KIC 8626021 by con- 
sidering six parameters of their DB white dwarf models: T^s, 
M», two parameters describing the C-O core composition pro- 
files (XcQ'fmfl and two parameters that define the envelope 
structure (Mgnv, A^Hefl They fix Mgnv and vary the remainder 
5 parameters. In contrast, in our models we have only 2 free 
parameters: T^s and M», and the chemical structure at the core 
and envelope is kept fixed according to the predictions of the 
evolutionary computations. In order to compare the quality of 
our best fit with the results of BK011, we compute the Bayes 
Information Criterion (BIC; Koen & Laney 2000): 

BIC^A^pJ^J + logcr^, (5) 

where A^p is the number of free parameters, and the number 
of observed periods. In our case, - 2 (stellar mass and ef- 

BK01 1 define Xq as the central O abundance and qt^^ as the edge 
of the homogeneous C-O core. 

^ BK011 define M^„y as the mass coordinate of the base of the He 
layer, and Mne as the mass coordinate of the point where the He abun- 
dance raises to 1 . 



of the DBV star KIC 8626021 



0.4 - 




30000 29000 28000 27000 26000 25000 24000 23000 22000 21000 



Fig. 5. The inverse of the quality function of the period fit in 
terms of the effective temperature. The vertical gray strip indi- 
cates the spectroscopic Tetf and its uncertainties (Teff = 24 900 + 
750 K). The curves have been arbitrarily shifted upward (with a 
step of 0.02). The strong maximum in ix^T^ corresponds to the 
best-fit model (model 2 in Table [B, with M, - 0.664 and 
Teff = 27 263 K. 

Table 2. Comparison between the observed periods of KIC 
8626021 and the theoretical periods of the best-fit model (model 
2 in Table[B. 



Hobs,/ 


A, 


H* 


e 


k 


Si 




[s] 


[ mma] 


[s] 






[s] 


[10-"* s/s] 


197.11 


2.92 


200.43 


1 


3 


-3.32 


5.60 


232.02 


5.21 


233.69 


1 


4 


-1.67 


8.66 


271.60 


1.87 


271.22 


1 


5 


0.38 


6.84 


303.56 


1.29 


301.42 


1 


6 


2.14 


9.91 






339.69 


1 


7 




6.67 


376.10 


1.05 


376.50 


1 


8 


-0.40 


12.13 



fective temperature). The smaller the value of BIC, the better 
the quality of the fit. We obtain BIC - 0.85, which is substan- 
tially larger than the BIC value of the best fit model of BK01 1 
(BIC = -0.41). This means that our period fit is somewhat 
poorer than theirs. Notwithstanding this, our asteroseismologi- 
cal model still provides a very satisfactory fit to the periods of 
KIC 8626021. 

The last column in Table [2] shows the rate of period change 
of the fitted pulsation modes. Our calculations predict all of the 
pulsation periods to increase with time (fit > 0), in accordance 
with the decrease of the Brunt-Vaisala frequency in the core of 
the model induced by cooling. Note that at the effective tempera- 
ture of KIC 8626021, cooling has the largest effect on n,t, while 
gravitational contraction, which should result in a decrease of 
periods with time, becomes negligible and no longer affects the 
pulsation periods. Until now, no measurement of li in a DBV 
has been assessed, although important efforts to measure the 
rate of period change in at least one star at the blue edge (EC 
20058-5234; Dalessio et al. 2010) are being currently carried 
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Table 3. The main characteristics of KIC 8626021. The second 
column corresponds to spectroscopic results (0EA11), whereas 
the third and fourth columns present results from the asteroseis- 
mological study of BK011 and from the best-fit asteroseismo- 
logical model of this work, respectively. 



Quantity 


0EA11 


BK011 


This work 


Teff [K] 


24 950 ±750 


29 200 


27 263 


M. [Mo] 


0.56 ± 0.03 


0.570 


0.664 


logg [cm/s^] 


7.91 ± 0.07 




8.099 


log(L./Lo) 






-1.14 


log(i?.//?o) 






-1.93 


A^Hctotal [Mq] 






3.6 X 10-3 






-2.80 


-1.63 


MHe <*> 




-6.30 


-5.95 


Xo <"> 




0.60 - 0.65 


0.65 


qim 




0.36 


0.46 



Notes: <*' Mass coordinate, defined as log(l - (see BK011). 

**' Central abundance by mass. Coordinate in units of M, (see 
BK011). 



out, and a possible determination of the rate of period change 
for other DBV star (PG 1351+489) has been reported (RedaelH 
et al. 2011). 

The main features of our best-fit model are summarized in 
Table [3] where we also include the parameters of KIC 8626021 
extracted from 0EA1 1 and the seismological model derived by 
BK011. In the Table, the quantity Mee.totai corresponds to the 
total content of He of the envelope of the model. In order to 
make the comparison easy, we include in the table the four pa- 
rameters employed by BK01 1 that define the chemical profiles 
at the core and envelope of their models. Note that BK011 do 
not specify the value of MHe.totai of their best-fit model. The loca- 
tion of the best-fit model for KIC 8626021 both according to our 
study and that of BK01 1 in the log Teff - log g plane is shown 
in Fig. [1] Both asteroseismic studies lead to the conclusion that 
KIC 8626021 should be closer to the blue edge of the DBV in- 
stability strip than spectroscopy suggests. 

5. Discussion 

The results of our asteroseismological analysis point to a higher 
stellar mass of KIC 862602 1 than predicted by spectroscopy. We 
arrive at such conclusion through both the period spacing and the 
individual periods exhibited by the star. Regarding the effective 
temperature, our work indicates a higher Teff than the spectro- 
scopic measurement. 

Our results also differ somewhat from those of the seismo- 
logical analysis of BK01 1 . Specifically, we obtain an asteroseis- 
mological model that is more massive and cooler than that of 
BK01 1 . The fact that in both independent analysis a good match 
to the observed periods is found can be understood on the basis 
of the asymptotic behavior of ^-mode pulsations, that predicts 
that a lower effective temperature is compensated by a higher 
mass. What is interesting and exciting is that, in spite of the sub- 
stantial differences in the white dwarf modeling (in particular, 
the quite different composition profiles, that lead to significant 
differences in the pulsation periods) both analysis agree that KIC 
8626021 is a hot DBV. This is in agreement with analyses based 
on the average period spacing and also with the fact that low pe- 
riod modes are present in KIC 862602 1's pulsation spectrum, as 
is also observed for the hot DBV EC20058. 



6. Conclusions 

In this paper we have presented a detailed asteroseismic analysis 
of KIC 8626021, the first pulsating DB white dwarf star discov- 
ered by the Kepler Mission, on the basis of the full evolution- 
ary DB white-dwarf models presented in Althaus et al. (2009a) 
which were computed for a wide range of stellar masses and He 
envelopes. These DB white dwarf models are characterized by 
consistent chemical profiles for both the core and the envelope. 
These chemical profiles are the result of the computation of the 
full and complete evolution of the progenitor stars from the zero 
age main sequence, including the core H and He burning phases, 
the thermally pulsing asymptotic giant branch phase, the born- 
again episode that is responsible for the H deficiency, and from 
time-dependent element diffusion predictions during the white- 
dwarf stage. 

By considering the mean period spacing exhibited by the 
star, we found that KIC 8626021 should have a stellar mass in the 
range 0.60 < M^/Mq < 0.87, substantially larger than those de- 
rived by previous spectroscopic (M* ~ O.56M0; 0EA1 1) and as- 
teroseismic (M, ~ O.57M0; BK01 1) studies. We also found that 
period-to-period fits point to an asteroseismological model with 
an effective temperature of ~ 27 300 K, in strong conflict with 
the spectroscopic estimate (Teff ~ 24900 K). Our results are in 
agreement with the recent asteroseismic analysis of BK011 on 
KIC 8626021, in particular regarding its effective temperature. 
In fact, these authors conclude that KIC 8626021 must be a hot 
DBV with Teff ~ 29 200 K. It would be interesting to see what 
a spectroscopic analysis based on higher signal-to-noise spectra 
will tell about the surface gravity and effective temperature of 
the star. 

If KIC 8626021 is a hot DBV, as first found by BK011 
and confirmed now by our results, then further monitoring of 
KIC 8626021 with Kepler in the next years probably will al- 
low a measurement of D, which in turn could open the possi- 
bility to constrain the plasmon neutrino emission rate (Winget 
et al. 2004; BK01 1). This endeavour is evidently dependent on 
the models, which can always be improved. Uncertainties in the 
models can also be assessed (Bischoff-Kim et. al. 2008). 
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